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The evolution of non-adiabatic perturbations in models with multiple coupled perfect fluids with 
non-adiabatic sound speed is considered. Instead of splitting the entropy perturbation into relative 
and intrinsic parts, we introduce a set of symmetric quantities, which also govern the non-adiabatic 
pressure perturbation in models with energy transfer. We write the gauge invariant equations for 
the variables that determine on a large scale the non-adiabatic pressure perturbation and the rate 
of changes of the comoving curvature perturbation. The analysis of evolution of the non-adiabatic 
pressure perturbation has been made for several particular models. 
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Multi-component models play an important role in modern cosmology. They are required to study a large scale 
structure formation [H, 0| and are used in multi- field inflationary models Multi-fluid models are also essential in 

the context of preheating [1, 0] • It is well known that in cosmological models with multiple fractions the growing mode 
' solution of the curvature perturbation on comoving hypersurfaces 1Z remains constant on a large scale in the absence 
of any entropy perturbations P, 0> 9 • This quantity is valuable to cosmological applications, because, for example, 
it allows one to relate the perturbations generated at a stage of inflation with the primordial perturbations in the 
$H \ early radiation dominated era. In the presence of the non-adiabatic pressure perturbations the comoving curvature 
W)i perturbation 1Z always changes with time. Thus, the non-adiabatic perturbations may affect the observable power 
spectrum of the adiabatic perturbations. 

A scalar field and a cosmological fluid seem very different, but the distinctions are not fundamental. Any barotropic 
^ ' perfect fluid can be represented as a scalar field with a nontrivial kinetic part @. On the other hand, in the linear 
0^ , perturbation theory the scalar field with a self-coupling potential can be treated as a perfect fluid with sound speed 
• c 2 s = SP/Sp (computed in the fluid rest frame) different from adiabatic sound speed c^ ad ^ [Tfj| . For this reason, an 
effective fluid description is often used. Under this approach, a scalar field can be fully described by specifying its 
sound speed and the equation of state, and by phenomenologically introducing an energy-momentum transfer between 
the scalar field and other fluids. In slightly different and more general formalism, any system of N coupled canonical 
scalar fields can be modelled by N "kinetic" fluids with a stiff equation of state interacting with one "potential" fluid 
I | with a vacuum equation of state fll| . 
• • . A valuable class of multi-component models is formed by the systems containing the coupled scalar fields and the 
barotropic fluids. These are, among others, the cosmological models with coupled cold dark matter and dark energy, 
considered as a canonical scalar field usually called quintessence. For these models there has been found the possibility 
of existence of large scale instabilities at the radiation dominated stage, caused by the fast growth of non-adiabatic 
perturbations |l2f. Such early time instabilities have received much attention in a variety of scenarios EMI. 

Although non-adiabatic perturbations were investigated long ago, there still remain some unresolved issues. Cur- 
rently, the definitions of the relative entropy perturbation for scalar fields (l6l . [l7| (note that the authors of Refs. 
[H, [H used adiabatic and entropy field combinations instead) and of general perfect fluids with intrinsic entropy 
perturbations (Tll.[2fj| are different. Thought Ref. [ll[ gives a unified definition of relative entropy that is appropriate 
both for barotropic fluids and for scalar fields, its use leads to some theoretical difficulties. In particular, even in the 
simplest single field models this method assumes, in general, a nonzero relative entropy between the "kinetic" and 
"potential" fluids, which seems somewhat artificial. In the phenomenological approach there are also some problems 
with the description of the interaction of these auxiliary fluids and others. It is also desirable to study the adiabatic 
condition in more detail. For example, in the recent paper (2l| an analytical attractor solution for dark energy per- 
turbations in the synchronous gauge at a constant dark ener gy e quation of state parameter was obtained. For this 



iquai 

solution, the usually imposed generalized adiabatic condition [18|, does not take place 
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The aim of this work is to, at least partially, fill the gaps. The paper is organized as follows. In Section [TT| we 
briefly review the perturbed Einstein and continuity equations for a general multi-fluid model with energy exchange. 
Assuming that the total energy density can be represented as a sum of energy densities of some perfect fluids, we 
write, in Section IIII1 an evolution equation for the comoving curvature perturbation 1Z and an expression for the 
non-adiabatic pressure perturbation P na d- Instead of using of relative entropy perturbations, we introduce a set of 
symmetric quantities Su that determine the P na d in models of considered type. In Section Mil we also derive the 
gauge invariant equations for variables that describe the evolution of non-adiabatic density and pressure perturbations. 
Some applications of the proposed formalism are discussed in Section IIVI In particular, the case of two minimally 
coupled fluids and applications to cosmological models with coupled dark matter and quintessence dark energy is 
studied in detail. We conclude the paper in Section IVl 



II. BASIC EQUATIONS. 

Let us consider scalar type perturbations at a FRW background. The general line element for scalar perturbations 

is 

ds 2 = a 2 (r) {-(1 + 24>)dT 2 + 2B^drdx l + [(1 - 2?A)% + 2E tij ]dx i dx i } , (1) 

where four scalar 0, ip, B, E are the first-order quantities and r is the conformal time. 

A perfect /-fluid with density pi — pi + Spi, pressure Pi = Pi + 5Pi, four-velocity Uj — dx^ /ds and vanishing 
anisotropic stress is defined to have an energy-momentum tensor of the form 

T£, = (p/ + PiK^ + P/<^, (2) 

The perturbed four- velocity can be written to the first-order terms as 

u i=~ (1 _ 0)> w j > u iu = a [-(1 + 4>), c^ + Bj], (3) 

Here, we introduce the velocity potential vi, since the fluid flow is irrotational for scalar perturbations. 

Throughout this paper, we will work in the Fourier space. The first-order perturbed Einstein equations yield 

3W. (if}' + H<fi) + k 2 ijj + k 2 H (E' — B) = -AnGa 2 Sp, (4) 

ij}' + U(j) = -4irGa 2 (p + P)(B + vi), (5) 

ip" + 2Hi// + He/)' + (2H' + H 2 )(t) = 4irGa 2 6P, (6) 

{E' -B)' + 2H(E' -B) -<t> + i> = 0, (7) 

where A: is a comoving wave number, the prime denotes the differentiation with respect to the conformal time t, and 
H = a' /a. 

In the general case of N coupled fluids divergence of the energy-momentum tensor gives 

r# = Q?, (8) 

The four-vectors Q^ A are related by the constraint [l[ 

JV 

EQ/=°< (9) 
i=i 

which results from the conservation law of the total energy-momentum tensor. For convenience, one can decompose 
these 4- vectors into two parts 

= Qiu» + Ff, Qi = Qi + SQi, u„F^ = Q. (10) 

Here is the overall four-velocity, Qi is the energy density transfer rate and Fj is the momentum density transfer 
rate of /-fluid in the total matter gauge. We have Fj = a _1 (0, fj 1 ) up to first-order, where /j is a momentum transfer 
potential [T3 |. 

The continuity equations for coupled perfect fluids can be obtained by the linearization of conservation equations 
([8|). As a result, one can write the background 

p'j = + wi)pi + aQi (11) 
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and the perturbed equations [12j 

V/ + 3H {8pi + SPj) - 3 (p 7 + Pi) ip' - k 2 (p z + Pj) («j + £') = + aSQi, (12) 

[(pi+^)(t>/+B)]'+4K(pj+P 7 )(t;/+B) + (p/+ft)0 + <yi^ = aQ I (v + B)+af I . (13) 

Pressure perturbations can be expressed in terms of density and velocity perturbations as (T3 | 

= C^Pl - i4l ~ C s/(ad))P/p- ( 14 ) 

where the velocity perturbation Oi = —k 2 (B + vi), c 2 Is ^ airj = Pj/p'j is the /-fluid adiabatic sound speed and c 2 j — 
SPi/Spj is defined in the /-fluid rest frame. The values of (? Is r a< j\ and c 2 gI are equal for the barotropic fluid (i.e. fluid 
for which the pressure depends only on the density), but this equality may not be satisfied in a general case. For 
example, for a canonical scalar field one has to set c 2 sI = 1 [lOfl . 
Using the notation 

qi = Qa ,,_ , _ ; (15) 
3H{pi + pi) 

and the identity 

w'j = (^\ = m(l + wi)(l - qi) (wi - c 2 sI(ad) ) , (16) 
the equations ()12 p .()13 p can be rewritten as 

(A, - 3V0' + (c 2 sI - c 2 sI{ad) + qi(c 2 sI{ad) + 1)) A, 

where A/ = 8pi/(pi + Pi). 

III. NON-ADIABATIC PERTURBATIONS. 



-9H 2 (c^ - 4 (ad) )(l - qi )tL + 6x- k 2 (E' -B) = m qi U + J , (17) 



The comoving curvature perturbation 1Z = ib — H(v + B) at spatially flat background can be represented in terms 
of the longitudinal gauge-invariant quantities [22] 

$ = 4> + H(B-E') + (B-E')', * =ip -H(B - E') (19) 



as 



2*' + %$ 
3 rt(l -I- w) 



The time derivative of equation (|20j) gives the well known equation 



2 1 f k 



.2 



ft'=o ?77i I „, A ^rAad-fc 2 c 2 ^ , (21) 



3 W (1 + w) { 2 " ' " aa " " s ( a 
where the non-adiabatic part of the pressure perturbation is defined by 

Pnad = SP- C 2 (ad) 5p. (22) 
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Substituting the expression c^ ad j = ^~ c 1i( a d) ^ or ^ ne overa ll adiabatic sound speed, at the arbitrary number of 
coupled perfect fluids, one can obtain 



Pnad = ^^{fsI- c2 sI{ad))^n+^{pi+Pl)c 2 sI qi^ 
I ? 



\- (pi + Pl)(pj + Pj) ( 2 iql7 / 2 2 \ n 

J^p \^sAi.i + (c s/ - c sJ(od) J (1 - qi)-p) , (23) 



i, J 

where A = Sp/ (p + P), A/j = A/ — Aj, and the gauge-invariant quantity 



e m = 8p + m{p + P)^ (24) 



is a comoving density perturbation [23]. 

The equation J23J) contains quantities A/j that are not gauge invariant at non-minimal coupling. Hence, following 
Refs. [11], [20], [24], it is convenient to use the gauge invariant variables 



6 pi bp j 



(l-qi)(pi + Pi) (1 - qj){pj + Pj) 



(25) 



(l-gi)(l-«j)A 7J (26) 



instead. The latter quantities are well defined even at qj = 1. 
In these new variables, the expression ()23|) can be rewritten as 



p srP'ifz 2 \ , 1 \ - (Pi + Pi)(pj + Pj) g (07 , 

Pnad = 2^J [ C sl - Csl(ad) J £ ™ + 3 <L J^p SlJ > ( 27 ) 

where 

S U = (4 - Ajj + ZH ((<& - C ^ (Qd) ) (1 - g/ ) - (f sJ - c 2 sJ{ad) ) (1 - qj )) (28) 
The perturbed Einstein equations Q and ((5]) lead to the relation 

^ = -4irGe m , (29) 

and in the long- wavelength limit k — > we obtain e m — > 0. Hence, the terms with e m in the equation (|27[) can be 
neglected on a large scale. 

The symmetrical quantities Su are closely related to the entropy perturbations and adiabatic condition. Using the 
definition of the scalar field energy-momentum tensor, the adiabaticity condition for systems of canonical scalar fields 

^-^ = (30) 

can be rewritten as 

% = 0, (3D 

where I, J = 1, N. 

The adiabatic mode in multi-fluid models with barotropic coupled fluids is determined by the conditions [lij 



0. (32) 



In both cases, quantities Su vanish. The case of systems with barotropic fluids and scalar fields is more complicated, 
since the leading contributions to Su depend on the gauge choice. In cosmological models with minimally coupled 
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quintessence dark energy, one can put Sjj — for the growing adiabatic mode in the longitudinal gauge, but some of 
quantities Sjj are different from zero at leading order in kr at "generalized initial adiabatic conditions" of (2l| in the 
synchronous gauge. The "generalized initial adiabatic conditions" provide only that the inequality P na d *C Sp holds. 
In any case, to find P na d and the comoving curvature perturbation TZ, it is necessary to know the values of Sij. 

The equation (|27[) shows that the non-adiabatic pressure perturbation depend on A/j and Qjj through their linear 
combinations Sij. This fact allows to simplify the analysis of the evolution of the non-adiabatic pressure perturbation 
that is important for the study of possibility of the non-adiabatic early instabilities in the models with coupled fluids. 

The quantities Su are completely determined by the values of Ajj and 6jj. Following Refs. [24],[20],[llJ, we 
introduce the notations 

. _ 5Qi Q'j afi afj ~ _ q/ej qjij 

Qi iUQi pi{l + wj) pj(l+wj) 1-9/ I-?,/ 

Then for the system of N coupled fluids (taking into account p = ^ 7 pi) the equations (|17|). (|18p yield 



A /, \ (j^ w (4 - 4 M) ) - m qi (cl I(ad) + 1)) E y ( A « + 

+\ ( j-^- - 3^ (4 c 2 SJ(ad) ) - m qj (cl J{ad) + i^Y,y + A ^) 

+ 7}K {fal - 4(ad) + C lj - C ij(ad) + <H{ C ll(ad) + l ) + <U{ C lj{ad) + *■)) A " 

1 ( I'l l'j \ A , 9 u2 / 2 2 ,2 2 \^PK+Pk (0 IK Ojk\ 

= (i^7 T^j) tf « - 3H - + - C ' J ) ftp' (34) 

-p- - ^{c s i-c s j-qi(l + c sI ) + qj(l + c s j))}_^ _ + p + -p- j 

+^ (l - §4 " §£, + §«,(! + 4) + + 4)) ^ 

-J (4a - ?/) - 4(i - E % + M 

-\ (4(1 - ?/) + 4(1 - ?./)) A 7J = ((1 - - (1 - gj )^) - /„, (35) 

where I,J,K = 1, ...,N and 

# uc = + fc 2 l (36) 
H 

is the velocity perturbation on uniform curvature hypersurfaces. 

These equations agree with the equations of Ref. [24] on a spatially flat background at a vanishing anisotropic 
stress 1 . On a large scale, after substituting the explicit expressions for Eab and Jab, the equations (|3"4T) . (f3"5|) form 
a closed system of equations. 

The system of equations (|34|). (j35]) contains N(N — I) independent equations for the antisymmetric quantities Ajj 
and Ojj. Meanwhile, for the analysis of the evolution of the non-adiabatic pressure perturbation, it is suffices to know 
only the values of N(N — l)/2 symmetric quantities Sij. In many cases of practical importance, on a large scale, it 
is not need to solve the complete system of equations (f34| and (|35|) , and one can derive and solve the equations for 
Sij. In what follows, we consider some applications of the new variables. 



1 There is a typo in the common sign of the second line of Eq. (2.196) of Ref. 24] . 
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IV. SEVERAL MODELS. 
A. Two non-coupled fluids. 



Models with non-interacting fractions are often used in modern cosmology, especially in inflationary models. They 
are of interest as toy-models when considering the radiation dominated universe. Primarily, we are interested in 
studying the influence of non-adiabatic sound speed on rate of decay of non-adiabatic perturbations in cosmological 
models with quintessence. We would also like to clarify the general arguments [25j that the non-adiabatic perturbations 
evolve independently of the adiabatic ones within the effective fluid formalism applied in the paper. 



1. Non coupled barotropic fluid and scalar field. 



Consider a simple model with two minimally coupled perfect fluids, indicated by subscripts A and B, where the 
-B-component is a barotropic fluid. In this case, the equations ([Mil . ([55)1 are reduced to 



A' AB +3H (c 2 A -c 2 



9' 

AB 



ft- 



sA L sA(ad) 
AB 



MB 



'AB 



k 2 k 2 

The definition ([28]) gives 



2 PB i „2 PA 



'sA -i 
P 



sB -, 
P 



MB 



3ft 



ve AB 

Iab 

k 2 



-3ft (c 2 , -c 2 



A(ad) 



p + P 



-sA ~ ^sB ) 



p + P 



Sab = {c 2 sA - c 2 sB ) A AB + 3ft (c 2 A - c 2 A{ad 



k 2 



(37) 
(38) 

(39) 



Here we restrict ourselves only to the case of constant c 2 A , c 2 B , wa, wb- At constant parameters wa, wb the 
equation (Til?)) implies the relations 

c2 sA(ad) = W A, C 2 sB{ad) = W B . (40) 

The equations (|37|) . ([38| can be combined to form the second order equation 



Q 



AB 



n 



1 i Q | „2 PB_ , 2 P_A 



2 2 
C sB ~ C sA(ad) 



3c; 



B " 



+3ft 2 



c 2 sA ^ + c 2 Pa 



-sB" 



H 2 



,2 PB i 2 Pa 



n 2 m 



A(ad) £' 
2 



sB 



AB 



-sB 



C sA C sA(ad) A 

n 9 JAB 



P+P 1 



where 



sA 



ad ) \ ( r 2 „2 

I [CsA c sA(ad) 



" u sA v sB , 

For two barotropic fluids with the indices A and B, the equation (|4"Tj) is reduced to 

S' 



S'ab + H 



| _ •! I r 2 Pb , 2 Pa 

1 J \ c sA -, T c sB -i 



, 2 I 2 PB , 2 PA \ q U 2 ( 2 2 \ 2 6 



(41) 



(42) 



(43) 



The additional multiplier k 2 on the right hand side of (1431) ensures that in the large scale limit k -+ 0, the source term 
is negligibly small, and the adiabatic perturbations do not affect the entropy evolution. 

When the A-component is a canonical scalar field, the implicit form of the equation (|41|) looks rather complicated. 
We write it only for the case of radiation and subdominant scalar field. Then, on a large scale, we obtain 



S'ab + 3ft fl - c 2 s A(ad)) S'ab + 6ft 2 (l - c 2 A ( ad })SAB — -ft 2 ^5 - 3c^4( a(i )J - (? s A( a d) 



p + P' 



(44) 
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Although this equation is gauge invariant, according to the definition of Sab, one can drop the source only if 
max{\Aji\, \/S.b\,H\Q a/ k 2 \,%\6 b / k 2 \] ^> This condition is not satisfied, in general, in the synchronous gauge. 
Indeed, for the "generalized adiabatic initial conditions" of Ref. [2l| 

a - - -f«X% {kT)2 ' (45) 



c 



2 4- 6w A + 3c 



sA 

(fcr) 3 fc, (46) 



we get 



Sab = (1 - «u) ^f^^rf, = -\^ { krf, (47) 

3(7 — ow A ) 2 p + P 3 2 

where C is a constant. A direct verification shows that (|4"4")) holds for the "generalized adiabatic initial conditions" 
only if we do not neglect the source in this equation. Thus, in the synchronous gauge, the adiabatic perturbations 
give rise to non-adiabatic perturbations even on a large scale. 

It is interesting to compare equation (|43p and (|4"4")l at the radiation dominated stage in case the quantities Sab are 
large enough for the sources to be neglected. The equation (|4"3")l yields on a large scale 

S% + H (1 - 3w A ) S' AB = (fluid -fluid). (48) 

The equation (|4"4")) takes the form 

S'a B +H(3-3wa)S A b+6'H 2 {1-wa)Sab=0 (scalar field - fluid) . (49) 

It follows that, if the A-fluid is the minimally coupled scalar field, the damping force is stronger. In addition, when 
wa < 1, the coefficient at Sab is positive. Hence, at negative wa, the quantity Sab rapidly approaches the asymptotic 
solution. For this reason, in numerical codes like, CMBFAST and CAMB, the initial dark energy perturbations are 
set by default to zero. 



2. Two scalar fields. 



The case of two non-interacting scalar fields, and the more general case of interacting fields are studied carefully 
by different methods. Here, we look at the second order equation for Sab on ly- From the original equations (I34|) and 
(1351) . it is easy to write an equation for Oab- It has the form 

^ + *(l-3(cL M )f + 

+n 2 ^3 - 2— - 3 (c? 8 A(ad)^- + C l B {ad)^f^j = 3 ^ ( C sA(ad) ~ C sB(ad)) - + p- ( 50 ) 

At constant adiabatic sound speeds, this equation can be rewritten in terms of Sab as 
s'k, + n (l - 2^ -3 (oXfli + c^ M |i)) S' + k*S AB 

= -9H 2 (c 2 A(ad) - Cg B ( od )J p + P' 

In general, one has to take into account the source term even on a large scale, but usually the scalar fields are 
considered in the longitudinal (E = B = 0) gauge. The well known long-wavelength adiabatic solution in this gauge 
[ltl can be rewritten as 

<t> = <^ = c(\-\ [ a 2 dr\ , (52) 



A A 


= A B 


e A 


0B 


k 2 


k 2 



Tl 



~n2 



a 2 dr, (53) 
a 2 dr, (54) 



Tl 



where C and t\ are constants. By the definition (|24[) . the value of e m vanishes now in the leading order, and the right 
hand side of the equation ([51]) can be neglected. 



B. Special case of three fluids. 

Consider the Universe filled with quintessence dark energy (A), dark matter (B) and radiation (C), and assume 
that the dark energy and the dark matter couple. Since the variables A/j, 9u are antisymmetric and constrained by 

A AB + A BC + A c a = 0, 9 AB + 8 BC + Oca = 0, (55) 

there are only four independent variables among them. 

Using the relations (|55|) and the equalities c 2 sB = c 2 sB ^ ad y c 2 c — c^ c , ad s, which are valid for barotropic fluids, the 
equations (f34|) and (|35|) can be reduced to set 



A' AC - ™ (cL - cl A(ad) ) - m qA {cl A{ad) + 1)) (^y^A AC - ~^A BC ) 

' p B + P B + pc + Pc Oac Pb + Pb Qbc ' 



+9H 2 [ ci A - d 



s A sA(ad)) y - +p k2 - + p k2 

Pb + Pb + Pc + Pc 1 , PA + P A \ a q_A Pb + Pb Q 

Vac — , f, Ubc 



P + P 1-qA P + P J l-q A p + P 

1A a Qa , / 2 2 \ e r, 

l^f™ ~ 6H \ C ° A ~ C ^)) P' 



= 3UE AC - t^-0 uc - m (c 2 sA - c 2 sA(ad) ) (56) 



y AC 



nj (t o 2 \ ®AC on, ( 2 2 /-, , 2 \\ ( PB + Pb + PC + Pc 9 AC PB + Pb BC 

H(l- Zc sC ) -nr-ZU [c sA - c sC -q A {l + c sA )) 



k 2 ' v sU) k 2 '" VsA sC *«\->-.AJjy p + p k 2 p + p k 

((1 - q A y sA clc) {^±^A AC - ^A SC ) - c 2 sC A AC = ((1 - q A )c 2 A c 2 c ) ^ - f AC , (57) 

Pb+Pb , Pa + Pa+Pc + Pc 1 \n Bc _PA±^ T 3B_ eAc = m ^ Bc _^ e ^ (58) 



p + p p + P 1 - q B J p + P 1 - q B " 1 — qB 

1 BC , v A o 2 \ °BC ™ / 2 2 ft , 2 \\ f PA + Pa + PC + Pc &BC PA + ?A AC 

^( 1 - 3c s c)^2--3H(c sB -c sC - te (l + c sB )) — -_ -= — _. 



fe 2 V P "'-V^ ^ " \^ p + p p p + P k- 

- ((1 - to )c 2 B - clc) (^Aa bc - ^A AC ) - cl c A BC = ((1 - q B )c 2 sB clc) ~ - /bc(59) 
VP P J p + P 

At the radiation dominated stage the dark energy and dark matter are subdominant, — < 1, — < 1, which 
allows to simplify the equations. Furthermore, we assume that we can neglect terms with in the long-wavelength 
limit. As a result, on a large scale, we obtain the approximate equations 



A' AC + [c 2 sA - c 2 sA(ad) + q A {c 2 sA(ad) + 1)J A AC + 9H 2 (c 2 sA - c 2 sA(ad) )(l -QA)-*r - 3W(1 - q A )E AC> (60) 

3' 



2/2 2 N^C 

"-jf+H(l- 3c 2 sA + 3q A (l + c 2 sA )) ^ - c 2 sA A AC - Jac i (61) 

A^ c + 3^qb(^ bM) +l)A B c + e B c+qB0 uc = m(l-q B )E BC , (62) 

^+H(l-3 C 2 B +3 te (l + C 2 B ))^- C 2 s A BC = -f BC . (63) 



In the following we explicitly use the values c 2 A = 1, c 2 B = 0, c 2 c = 1/3. 
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1. Example 1. 



As a simple model, consider the coupling [2f 



Q 1 * = 1PaPb{u b - u A ), 



(64) 



where 7 is a constant. 

For such interaction, all the quantities Q A are zero and the evolution of background variables is the same as at 
minimally coupled fractions. Now q A = qs = = (B — and 



AC 



aiPB Vab 
1 + WA k 2 



BC 



-ai PA- 



JAB 



The equations (|62p and give 



A', 



BC °BC 



= o, 



9 'bc 
k 2 



n 



k 2 



0AB 
k 2 ' 



(65) 

(66) 
(67) 



We consider here only the case of wa < —1/3. Since for minimally coupled fractions pa oc a 3 ( 1+w a) and | a !£ A | <C 1 
at early times, we find that one can consistently assume that Abc — Obc — if Oab are not very large. Then the 
first two equations, (j6"0)) and (|61|). are simplified to 



A A c + 3H(l-WA)AAc + 9n 2 {l~w A )^§- = 0, 



9'ac 



n 



aiPB 



k 2 ' '"\H(l+w A 
These equations can be rewritten in the equivalent form 

ajp B 



- 2 



k 2 

A AC = 



s'a C +hii 



Sa 



c 



n 



H(l + w A ) 

A' AC + n(i-3w A )A AC = ~ms AC - 



5 2 ajp B 

Wa H ; r 

3 3H(l + w A ) 



A AC, 



(68) 
(69) 

(70) 
(71) 



Concerning the sufficiently early Universe, we have 
order equation 



H(1+wa) 



3> 1. Under this condition one can obtain the second 



a lPB x, Qa , 2 n \ a lPB 

-Sac + 3U {1-wa)- 



H(1+w a ) 



H(l + w A ) 



Sac = 0. 



(72) 



The presence or absence of instabilities is completely determined by the sign of the coupling constant 7. At negative 
7, there are an anti-damping force and tachyonic instability, which increase indefinitely in the limit a — > and lead 
to a catastrophic growth of Sac and a non-adiabatic pressure perturbation P na d- 



Another simple coupling is Qa 
densities are 



This background solution yields 



2. Example 2. 

-Qb = IPb At the radiation dominated stage the background dark sector 



PA 



:Pb, Pb oc a . 



3w A + 2 



(73) 



Eac ~ Qa^bc, Ebc ~ 0, 



(74) 



(75) 
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where 

Zw A + 2 
3(1 + w A ) 

The equations (f60"T) -(f63 |) reduce now to 

A' AC + 3n(l-w A + q A (w A + l))A AC + 9n 2 (l-w A )(l-q A )^- = 2,Uq A A BC , (77) 

— +n(6 qA -2) 1 ^-A AC = - TT ^ 1 ^ (78) 

and 

A' BC + 6 BC = 0, (79) 



k 2 > = °- ( 8 °) 

From the last two equations it follows that one can consistently assume that A B q = Obc = 0, if this condition was 
fulfilled initially. The first two equations (1771) and (178")) then give 

A' AC + m{l - w A + q A (w A + l)) A AC + 9H 2 (l - w A )(l - Qa)^ = 0, (81) 

9 -ML+H(6q A ~2) 9 -^ = A AC . (82) 

Using the definition of S A c, the equation (|76p . and substituting H = 1/r, one can write now the second order 
equation 

-„ llktu + 9-, 2 11^ + 9 5 

+ 7TT^a~ 5ac + ^TT^Z 5 ^ " °- (83) 

Assuming a power-law form of the solution S A c oc r™ , we obtain the algebraic equation for the power-law index 

t , 2+ 1Q Wa + 8, +2 11^ + 9 =0; 

1 + WA 1 + ^ 

with the roots 



^ = -^±^0. (85) 
1 + w A 1 + w A 

The equations (|2"7|) and (1731 imply that on a large scale P na d ~ ^-S A c oc t 3 Sac- Hence, P„ ac ; = Cit"+ + C2T 71 - , 
where Ci , C2 are constants and 



n ± = -^±i±VME!. ( 86) 

1 + 1 + W A 

The resulting expression agrees with the equation (73) of [12| and a further analysis is identical to the one given in 
Ref. EH. 



V. CONCLUSION. 



In this paper we considered the evolution of the non-adiabatic perturbations in the models with multiple interacting 
fluids with c 2 j ^ c s/( a ci)- We wr °t e the gauge invariant equations for the variables that determine the non-adiabatic 
pressure perturbation and the rate of changes of the comoving curvature perturbation 1Z in the models with energy 
exchange . The analysis of these equations was made for several particular models. One can see that the quantities 
Su, introduced in the paper, allow to clarify the features of the non-adiabatic perturbations evolution in multi-fluid 
cosmological models. 
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